Symmetries of multifractal spectra and field theories of Anderson localization 
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We uncover the field-theoretical origin of symmetry relations for multifractal spectra at Anderson 
transitions and at critical points of other disordered systems. We show that such relations follow 
from the conformal invariance of the critical theory, which implies their general character. We also 
demonstrate that for the Anderson localization problem the entire probability distribution for the 
local density of states possesses a symmetry arising from the invariance of correlation functions of 
the underlying nonlinear a model with respect to the Weyl group of the target space of the model. 
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More than half a century after its discovery, Ander- 
son localization [l| remains a vibrant research field. One 
of the central research directions is the physics of An- 
derson transitions 0, including metal-insulator transi- 
tions and transitions of quantum Hall type (i.e. between 
different phases of topological insulators). Apart from 
electronic conductors in semiconductor structures, exper- 
imental realizations include localization of light [3] , cold 
atoms [3] , ultrasound @ , and optically driven atomic sys- 
tems [6]. On the theory side, the field received a strong 
boost through the discovery of unconventional symme- 
try classes and the development of a complete symmetry 
classification of disordered systems 0,0-Q. These classes 
are characterized by additional particle-hole and/or chi- 
ral symmetries. Examples include disordered supercon- 
ductors and graphene. 

A remarkable property of Anderson transitions is the 
multifractality of wave functions, describing their strong 
fluctuations at criticality. Specifically, in d dimensions, 
the wave function moments show anomalous multifractal 
(MF) scaling with respect to the system size L, 



L d (\yb(r)\ 2q ) (xL- 



r q = d(q - 1) + A, 



(1) 



where (...) denotes the disorder average and A 9 are 
anomalous MF exponents distinguishing the critical 
point from a simple metallic phase (where A q = 0). 
Closely related is the scaling of moments of the local 
density of states (LDOS) p(r), 



(P q ) oc L~ 



= A q +qx p , 



(2) 



where x p = xi controls the scaling of the average LDOS, 
(p) oc L~ Xp . First steps towards experimental determi- 
nation of MF spectra have been made recently [!, [l(| ■ 
In Ref. [H}, an exact symmetry for MF exponents 



A„ = Ai_ 



(3) 



was derived for any critical system in the conventional 
Wigner-Dyson (WD) classes as a consequence of a more 



general relation 12, 13| for the LDOS distribution func- 
tion (and thus, for the LDOS moments), 



v(p) = P - a v(p- L ), ( P q ) = (p 1 -*) 



(4) 



Equation (j4]) is exact at the level of the nonlinear a model 
and is fully general otherwise; i.e., it is equally applica- 
ble to metallic, localized, and critical systems. While in 
general a models are approximations to particular micro- 
scopic systems, Eq. Q is exact in view of universality of 
the critical behavior [llj] . See also [14| . 

The goal of the present work is to reveal the field- 
theoretic basis underlying the symmetry relations (J3J and 
(|4]), and to generalize them to a broader class of systems. 
First, using arguments arising from conformal invariance 
at the transition, we show that relations analogous to 
([3]) are valid for a wide class of critical points in disor- 
dered systems (that need not be Anderson transitions) 
characterized by multifractality. Second, focusing on the 
problem of Anderson localization, we demonstrate that 
Eqs. ([3]) and (|4|) are manifestations of the Weyl group 
symmetry of the nonlinear cr-model theory. Finally, we 
use this to generalize Eqs. (j3]) and (j4]) to the unconven- 
tional symmetry classes C and CI (in the notation of 
Refs. [7, 8]). At the end we describe applications of our 
results to a number of specific disordered systems. 

We begin by presenting a general argument based on 
(global) conformal invariance of a critical system in d 
dimensions. Consider a system at criticality character- 
ized by operators O q representing moments of an observ- 
able of interest. In the case of Anderson localization on 
which we focus, this observable is the LDOS, and O q 
corresponds to p q , but one can apply the argument to 
a broader class of systems. Generically, the spectrum of 
the scaling dimensions x q of the operators O q is convex, 
x" q < (primes denote derivatives with respect to q) , sat- 
isfies xo — 0, and becomes negative at sufficiently large 
(positive or negative) values of q 0. Therefore, there is 
a single point g» (in addition to q — 0) such that x q ^ = 0. 
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Let us show that g* > 0. Indeed, it is easy to see that 
the derivative x' = {dx q /dq) q= o = ao — d + x p , where 
a Q = (dr /dq) g —o = Tq controls the scaling of a typical 



wave function amplitude, 



typ 



L a ° . Normaliza- 



tion of the wave function implies that ao > d Q. In 
the WD symmetry classes where x p — 0, this guarantees 
that x' > 0. In the unconventional classes, x p may be 
nonzero, with either sign. However, generalizing the con- 
formal invariance argument from Ref. [Hf, we can show 
that aQ — d+x p determines the typical localization length 
in a quasi- ID geometry, implying again that x' > 0. It 
follows immediately that q* > 0. 

According to the definition of the operators O q , their 
operator product expansion has the form [16| 



O p {n)O g {r 2 ) ~ |n-r 2 | 



c "0 



p+q 



In general, the operator O p + q has nontrivial scaling with 
the system size, (O p+q ) oc L^ x "+ q ; the existence of neg- 
ative scaling dimensions distinguishes disordered critical 
points from conventional "unitary" conformal field the- 
ories. However, in the case of p = q* — q, the oper- 
ator on the right-hand side of ([5]) has zero scaling di- 
mension x qr = 0. Therefore, the correlation function 
(Oq,-q{ri)O q (r2)) docs in fact not depend on the sys- 
tem size L (i.e. on the infrared regularization of the the- 
ory) . This allows us to apply the standard argument from 
conformal invariance 17] according to which the nonva- 
nishing two-point correlation function appearing in the 
expectation value of ([5]) implies that the dimensions of 
the nonderivative (T^j operators O qm - q and O q are equal; 
i.e. 

Xq X q ^_ q . (6) 

This is the generalized symmetry relation for MF expo- 
nents. Note that in the case of Anderson transitions the 
symmetry holds in general for the LDOS exponents x q 
(the scaling dimensions of local field operators) rather 
than for the wave function exponents A q . (For the WD 
classes x p — 0, so that x q = A q .) The symmetry point 
(?*/2 remains unspecified by the above argument. 

Below we show that a stronger symmetry relation be- 
tween moments of the LDOS, valid for nonlinear a models 
describing disordered systems of non- interacting fermions 
in the WD classes (A, AI, All) as well as those in the 
Bogoliubov-de Gennes classes with preserved spin rota- 
tion invariance (C, CI), 



(P q ) = {P q '- q ), 



(7) 



has a group-theoretic origin. As a result, the symmetry 
point g*/2 for these systems is determined solely by the 
symmetry class and is independent of further details of 
the problem (e.g., spatial dimensionality, presence or ab- 
sence of topological order, and whether the system is in a 
metallic, insulating, or critical phase). Before presenting 



the proof, let us first ask the following question: Assum- 
ing that we know that ([7]) holds with g* depending on the 
symmetry class only, what is a simple way to find q*7 It 
turns out that it suffices to analyze a zero-dimensional a 
model, equivalent to a random matrix (RM) model. 

To see this, let us consider a RM model and introduce 
a small broadening i <C A (where A is the mean level 
spacing) for all levels. A level with energy e gives a con- 
tribution ~ S/(e 2 + S 2 ) to the LDOS p at zero energy. It 
is easy to see that for small S the zero-energy LDOS will 
be governed by the level closest to zero. Thus, we get 



{pi) cc / de[S/(e 2 + 6 2 )]«P(e), 



(8) 



where P(e) is the distribution of the lowest energy level. 
For RM ensembles one has P(e) oc \e\ mi , where mi is the 
multiplicity of the long roots for the symmetric space of 
Hamiltonians. For all WD ensembles mi = 0, for class 
C mi = 2 and for class CI m; = 1. Equation ([8]) yields 
{pi) cx 5 m i +1 -4 for q > (mi + l)/2 and {pi) oc 5* for 
q < {mi + l)/2. This fixes the value of g* in the relation 
0: 



q* = mi + 1 = 



WD classes, 
class CI, 
class C. 



(9) 



We turn now to the derivation of Eq. . The a mo- 
dels are defined on symmetric superspaces G/K, where 
G is a Lie supergroup and if is a compact subgroup fixed 
by a Cartan involution 0, H, Ell • We focus first on the 
unitary WD class (A); in this case the a model field 
Q(r) = ,9(?')o'3ff(^) _1 is a 4 x 4 supermatrix satisfying 
Q 2 = 1. Here g(r) £ G and (73 is the third Pauli matrix 
in the retarded-advanced (RA) space. The moments of 
the LDOS at a point tq are given by [l2| 



{p")= DQ -(Q11-Q22+Q12-Q21) 



bb 



(10) 



where Q = Q(ro), with indices 1,2 referring to the RA 
decomposition and b,f to the boson-fermion one. The 
factor e - - 7 ^) results from integrating out Q{r) at the 
points r ^ tq and generically breaks the symmetry from 
G to K as a result of coupling the system to metallic 
reservoir(s). The only important property of the func- 
tion J-(Q) is its invariance with respect to the group K, 
i.e. F(kQk- 1 ) = F(Q) for any k £ K. This follows 
from the corresponding invariance of the action of the a 
model, including the boundary terms oc Str a^Q{r) ap- 
pearing at points r coupled to leads. In Ref. [12j the 
LDOS distribution function corresponding to (fT0| was 
evaluated by using the "standard" (introduced by Efe- 
tov) parametrization of the Q field [19], which led to 
Eq. (j4j. In order to uncover the group-theoretic basis of 
the symmetry, we will use an alternative parametrization. 
It is based on the Iwasawa decomposition for symmetric 
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superspaces 2(| 21 1 generalizing the corresponding con- 
struction for classical noncompact symmetric spaces [22| . 
By this decomposition, every element g £ G is repre- 
sented as g = nak with n £ N, a £ A, and A; £ K, where 
A is a maximal Abelian subgroup for G/K and N is a 
nilpotent group. The decomposition is unique, once the 
set of positive roots is fixed. (The corresponding root 
vectors form the basis of the Lie algebra of TV.) 

It is convenient to switch to Q — Qa^ and perform 
a unitary rotation Q — > Q = UQU^ 1 by the matrix 
U = (1 + i<r% + i(T2 + ias) /2 in the RA space, which cycli- 
cally permutes Pauli matrices: UujU^ 1 = Gj-x- The 
combination of Qij entering Eq. (|10[) then becomes 

(1/2)(Q U - Q 22 + Qi 2 - Q 21 ) bb = Q 22 , bb . (11) 

The Iwasawa decomposition of g leads to Q = 
na^u^n" 1 ^^, where we used ka^k^ 1 — 03 and aa^a^ 1 = 
a 2 a^. Upon the rotation Q — > Q, this takes the form 
Q = fWkG<xfiT (72, or explicitly 



Q = 
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(12) 

where * denote some nonzero matrix elements of nilpo- 
tent matrices. The variables x and y (which corre- 
spond to Ai = cosh 0i and A 2 = sin6* 2 in the standard 
parametrization fl9j ] ) parametrize the Abelian group A. 
This group is noncompact in the x direction and compact 
in the y direction. It follows from (TT2"j) that the matrix 
element (|TTj) is equal to Q 22 ,bb = e~ 2x . The integral (fTU)) 
for the LDOS moments thus becomes 



(P q ) 



DnDa e 



-2p(lna) — 2qx — Tina a 3 n <r 3 ) 



(13) 



NA 



where Dn and Da = dxdy are the invariant (Haar) mea- 
sures on N and A, respectively. The factor e - 2 P( lna ) is 
the super- Jacobian, with p(ln a) being the half sum of 
positive roots; for the present case 



P : 



ly. 



(14) 



Next, we perform the n integration involving only the 
last factor in the integrand of (fTB"]) . For this purpose, we 
use the Harish-Chandra integral theorem stating that for 
a if-invariant function [f(g) — flkgk" 1 ) for any k £ K] 



N 



dnf(na) = e p( - lna h f (a), I f {a w ) = I f {a). (15) 



The central point is that the function If {a) is invariant 
with respect to the action w : a — > a w by any element 
w £ W of the Weyl group W of G/K. The classical 
version of the theorem (|T5j) can be found, e.g., in (23[; the 
supersymmetric generalization (that we actually need) 



has been developed very recently [21] . The Weyl group 
acts on the Lie algebra of A; its elements are reflections 
with respect to hyperplanes orthogonal to roots. The 
element that will be important for us here is the reflection 
x ->■ -x. Substituting Eqs. (O and ([15]) into ([131), we 
obtain 



(jfl) = / dxdye^ 2 ^- i yi T {x,y). 



(16) 



Finally, by using the symmetry of Ijr^x^y) with respect 
to x — > —x, we obtain Eq. (0|. 

Thus, the symmetry relation (j3[) previously derived for 
the WD classes is a consequence of the Weyl group invari- 
ance. We can now extend this result to two new classes, 
namely, C and CI. By inspecting the above derivation, 
one can see that the value of <?* = 1 was determined by 
the coefficient in front of x in the half sum of positive 
roots p, Eq. (fj"4")) . The analogous formulas for the classes 
C and CI read p = — 3x + 2iy\ + iy 2 and p = —2x + 2iy, 
reproducing the values of obtained from the RM ar- 
gument, Eq. (JSJ). 

What about the remaining five symmetry classes? The 
above derivation based on the Weyl group invariance is 
not directly applicable to them because of a more compli- 
cated structure of the a model target space. Specifically, 
that space contains an additional U(l) factor in the case 
of the chiral classes AIII, BDI, and CII, and a 0(1) = Z 2 
factor for the classes D and DHL More work is needed 
to explore the peculiar physics of Anderson localization 
in these symmetry classes. 

The obtained symmetry relations are confirmed by a 
large body of analytical and numerical results for vari- 
ous disordered systems. For the WD classes, supporting 
evidence based on 2 + e expansion and simulations of a 
power-law random banded matrix model was presented 
ll| . Since then, a wealth of numerical results on 



in Ref. 

the 2D Anderson transition in the symplectic class All 
[24L the integer quantum Hall transition (unitary class 
A) [25], as well as the 3D transition in the orthogonal 
class AI (2(| have corroborated the relation ([3]). 

A thoroughly investigated representative of symmetry 
class C is the 2D spin quantum Hall (SQH) transition. 
For this system, it was proven analytically that, in the 
bulk, x p = 1/4, A 2 = -1/4, A 3 = -3/4 [13, In 
combination with the trivial values Ao = Ai = this 
yields xq = X3 = 0, and x\ = £2 = 1/4, in agreement 
with the symmetry relation ([6]) with g>* = 3. Further- 
more, these exponents have also been found for the SQH 
surface multifractality, with the results xo — £3 = and 
xi = X2 = 1/3 again respecting the symmetry. Fi- 
nally, the MF spectrum at the 2D SQH transition was 
studied numerically, both in the bulk and at the surface 
fjjl l30l ]. When expressed in terms of x q , the data per- 
fectly agree with the symmetry relation (JH]). 

There also exists a model in class CI that has been 
studied in detail. This is the model of 2D Dirac fermions 
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coupled to a random SU(2) gauge potential which is de- 
scribed by a Wess-Zumino-Witten theory [3l| and repre- 
sents physics at the surface of a disordered 3D topological 
superconductor 32[. Critical exponents for this model 
are known exactly: x p — 1/4 and r q = 2{q — 1)(1 — q/8) 
33], so that the LDOS MF spectrum reads x q = q(2 — 
q)/4, which clearly satisfies Eq. ([6]) with g* = 2 [in agree- 
ment with ©]. 

Further support for our results is provided by the 
LDOS distribution of a quasi-lD system in contact with 
a metallic reservoir. Far away from the contact the dis- 
tribution of In p is known to be Gaussian, with the ratio 
var(ln p)/{— hi p) equal to 2 for the WD classes, 2/3 for 
class C, and 1 for class CI [13, l35| . By calculating the 
moments (p q ), we recover Eqs. (0 [9|) for these classes. 

For the remaining five classes of Anderson localization, 
as well as for critical MF systems of other origin, our 
general arguments based on conformal symmetry predict 
a weaker (valid at criticality only) relation ([6]) to hold, 
with having the same degree of universality as criti- 
cal exponents normally have (i.e. they are controlled by 
a particular fixed point rather than solely by the sym- 
metry class). Several exactly solvable problems confirm 
this. In particular, the model of 2D Dirac fermions in an 
Abelian random vector potential (residing in chiral class 
AIII) has j3(| a parabolic MF spectrum with the sym- 
metry point q*/2 = 1/(4^) depending on the disorder 
strength qa- The same applies to the non- Abelian SU(A) 
(with N ^ 2) generalization of this model [33| (belong- 
ing to class AIII as well) where the symmetry point is 
N dependent: q^/2 — N/(2N — 2). Another example is 
the MF spectrum of the Ising disorder variable at the 2D 
Nishimori critical point, for which g» = 1 37]. Further- 
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more, the MF spectra of the harmonic measure of critical 
can be obtained by introducing (conformal) 



curves 

primary operators characterized by charges a (analogous 
to our q) [39(; their conformal weights h a (analogs of our 
x q ) form a parabolic spectrum with a symmetry point a* 
depending on the central charge of the model. 

To summarize, conformal invariance arguments reveal 
the general character of the symmetry relation (|6|) for 
critical disordered systems with MF scaling of moments 
of observables (represented by local field operators O q ). 
For the Anderson localization problem in the WD, C, and 
CI classes a stronger symmetry relation (J7]) holds, which 
is based on the Weyl group symmetry of the a model tar- 
get space. This stronger relation is not restricted to criti- 
cality, and is characterized by a symmetry point g*/2 de- 
pending on the symmetry class only. Future work should 
clarify the impact of the Weyl group symmetry on scal- 
ing dimensions of other composite operators, and explore 
the role of the U(l) and 0(1) degrees of freedom in the 
chiral, D, and Dili classes. 
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